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ABSTRACT 



in 

(N 
1^ 

^' 

^^ We present a study of the star formation (SF) region G38.9-0.4 using publicly avail- 

1-H able multiwavelength Galactic Plane surveys from ground- and space-based observa- 

J> tories. This region is composed of four bright mid-IR bubbles and numerous infrared 

k>( dark clouds. Two bubbles, N 74 and N 75, each host a star cluster anchored by a 

?-H single 09.5V star. We identified 162 young stellar objects (YSOs) and classify 54 as 

stage I, 7 stage II, 6 stage III, and 32 ambiguous. We do not detect the classical sign- 
posts of triggered SF, i.e., star-forming pillars or YSOs embedded within bubble rims. 
We conclude that feedback-triggered SF has not occurred in G38.9-0.4. The YSOs are 
preferentially coincident with infrared dark clouds. This leads to a strong correlation 
between areal YSO mass surface density and gas mass surface density with a power law 
slope near 1.3, which closely matches the Schmidt-Kennicutt Law. The correlation is 
similar inside and outside the bubbles and may mean that the SF efficiency is neither 
enhanced nor supressed in regions potentially influenced by stellar feedback. This sug- 
gests that gas density, regardless of how it is collected, is a more important driver of 



SF than stellar feedback. Larger studies should be able to quantify the fraction of all 
SF that is feedback-triggered by determining the fraction SF, feedback-compressed gas 
surrounding H II regions relative to that already present in molecular clouds. 

Subject headings: Stars: formation, Stars: pre-main sequence, Stars: protostars, (ISM:) 
HIT regions, ISM: bubbles 



1. Introduction 



Elmegreen & Lada (1977) proposed a mechanism where massive (>8 Mq) OB stars could 
trigger a new wave of star formation (SF) that would not have occurred otherwise. This 'collect- 
and-collapse' (CC) model showed that the expansion of an H II region into a uniform interstellar 
medium (ISM) can sweep up shells dense enough to collapse and form stars. A second proposed trig- 



gering mechanism is through radiation-driven implosion (RDI; Bertoldi 1989). Lefloch & Lazareff 



( 1994 ) demonstrated that ionizing ultraviolet (UV) radiation impinging upon a cloud can compress 



clumps of ISM to induce collapse. In their simulation, UV photons were directed at a preexisting, 
gravitationally stable core (n = 1550 cm^'^) embedded in neutral ISM [n = 10 cm~^). The re- 
sulting gas formed pillar-like structures with gravitationally unstable, SF cores at the tips, similar 



to what is observed in many star forming regions (SFRs) (Reach et al. 2004 Indebetouw et al. 



2007 Billot et al.||2010 ). Hydrodynamic simulations by Gritschneder et al. (2010) investigated the 
effect of radiation on molecular clouds with turbulent gas. Their "radiative round-up" (RR) model 
produced star forming pillars from the turbulent gas with distinct differences in line-of-sight gas 
profiles compared to RDI. 

The observational signatures of triggered SF should differ depending on the dominant mode, 
but in general, evidence for triggered SF will show patterns in the distribution of molecular gas, 
the distribution of YSOs, and the ages of YSOs relative to the H II regions. CC is a relatively 
slow process that takes on the order of a few Myr, which creates massive clumps that can collapse 



to form massive stars or star clusters (Deharveng et al. 2005). The expansion of the H II region 



triggers SF, which implies that YSOs should be found within, or adjacent to, the clumps that form 



on the periphery of the bubbles (Deharveng et al. 2008 Watson et al. 2008). The CC model also 



predicts that (in a uniform medium) the collected shell will fragment at regular intervals and lead to 
uniformly distributed SF clumps ( [Deharveng et ar]|2005 ) . As the ionization front pushes outward, 
new stars should continue to form (provided there is enough molecular material) and leave behind 
older, more evolved YSOs. Observationally, this predicts a YSO age gradient that can be used 
to identify cases of triggered SF. On the other hand, RDI can form young stellar objects (YSOs) 



in just a few tenths of a Myr, and the stars are intermediate mass, 1-7 Mq (Lefloch Sz Lazareff 



1994). RR differs from typical CC models which assume a static, uniform medium, and unlike RDI 



models, does not require pre-existing cores of molecular material. The RR model forms sub-solar 
mass protostars on timescales similar to RDI (~few 10^ yr). The stars formed in both CC and 



RR should have near-zero velocities relative to the expanding H II shell because they formed out 
of ISM already in motion. In RDI, the cores are formed before the expansion of the H II region, 
and therefore the resultant YSOs may have non-zero velocities with respect to the H II shell. The 
uniqueness of individual SFRs make it difficult to attribute specific features to a given mode of 
triggered SF (CC, RDI, or RR), however, the observational signatures mentioned above can be 
used to identify potential sites of triggered SF. 

This paper is the first in a series that will investigate triggered SF over a variety of Galactic 
environments. Our goal is to examine, in detail, the properties of the SFRs themselves to determine 
the input energy required for triggered SF to be an effective and efficient mechanism. This requires 
knowledge of the triggering and triggered populations, as well as the molecular cloud out of which 
they formed. To answer these questions we have gathered survey data from numerous sources, 
from the near-IR to the radio, to investigate the Galactic SFR G38.9-0.4. This region was chosen 
because it is relatively nearby (2.7 kpc) and lies along an unconfused line of sight. It also appears 
to host bubbles at different evolutionary stages and may provide insight into the growth of H II 
regions and the potential for triggered SF. 

Figure [l] shows a three-color image of G38.9-0.4 with UKIDSS K-band in blue, Spitzer IRAC 8 
/im in green, and Spitzer MIPS 24 fim in red. The green circles and text mark the main features of 
the region and indicate the apertures used for estimating the total gas mass and radio continuum 
flux associated with each feature in subsequent sections. This region contains two prominent IR 



bubbles (N 74 and N 75; Churchwell et al. 2006), two cataloged star clusters (#15 and #16, hereafter 



MCM 15 and MCM 16; [Mercer et al.|2005[ ), an intermediate-mass SFR (IMSFR, IRAS 19012-^0505; 



Arvidsson etrai:]|201()| >, an infrared dark cloud (IRDC) MSXDC G038.95-0.47 ( [Rathborne etaL 
2006, hereafter just IRDC), and an "ultracompact embedded cluster" (UCEC #8; Alexander &; 
Kobulnicky||2012 ). Magenta circles outline the areas used to find possible stellar clusters contained 
within the bubbles. The white crosses indicate the center positions of the features mentioned above 
as they are given in the cited references. Due to the IRAS beam size, the IMSFR is a blend between 
the mid-IR emission from MCM 15 and the N 75 bubble. The UCEC is the bright point-like feature 
located between the IRDC and N 74. 

In Section [2| we characterize the nature of the ISM and present evidence for the presence 
massive stars. Section |3] describes the method for finding YSOs. Section 4.1 assesses the main- 



sequence (MS) stellar content of the region. Section H^ is a breakdown of individual subregions. 
Section [5] is the discussion, including possible relationships between individual YSO masses and 
gas surface density, and between YSO mass surface density and gas mass surface density. We also 
discuss the YSO mass function in relation to other studies and the possibility of triggered SF. A 
summary of conclusions is presented in Section [6] 



2. Interstellar Environment 

2.1. Molecular Gas in ^^CO 

One indicator of active and future SF is the presence of molecular gas. To characterize this 
component we use data from the Boston University Five College Radio Astronomy Observatory 
Galactic Ring Survey (BU-FCRAO GRS; [Jackson et al.|2006[ ), which targeted the ^^CO J = 1 ^ 



transition at 110.2 GHz with a beam FWHM of 46" and velocity resolution of 0.21 km s~^. ^'^CO 
is much less abundant than H2 or ^^CO and thus remains optically thin at higher densities. We 
extracted spatially averaged velocity spectra from the GRS data cubes within N 74, N 75, cluster 
MCM 15, and the known IRDC. Figure [2] shows the ID velocity spectrum for each region with 
MCM 15 as the black short-dashed line, N 75 as the orange solid line, N 74 as the blue dot-dashed 
line, the IRDC as the green dotted line, and MCM 16 as the magenta, long-dashed line. The ^^CO 
brightness temperature peaks between 5.6 K and 10 K over the velocity range Vlsr = 39-42 km s^^ 
for all four regions, indicating abundant, probably related, molecular material over a small velocity 
interval. The peak velocities for each region increase slightly from north to south with a spread of 
about 3 km s~^ between the peaks. There is also low- level emission detected from 50 - 65 km s^^ 
for each region. There is no morphological similarity between this emission and the IR features so 
it is likely to be unaffiliated diffuse ISM. 

Figure [3] is a six-panel ^^CO channel map for the SFR from 35 to 45 km s~^, averaged every 
2 km s~^, overlaid on a GLIMPSE 8 /xm image. There is obvious spatial agreement between the 
molecular gas and IR emission. There is a strong peak coincident with MCM 15 at 36 - 38 km s~^ 
and a smooth transition across the image until another strong peak appears over the IRDC at 
43 km s~^. The 43 km s~^ panel also seems to show a ^^CO hole over the center of N 74. This 
feature may be evidence for a massive star sweeping out molecular material. The smooth variation 
in ^^CO velocity across this region suggests that all of these features are affiliated with the same 
molecular complex at a common distance. We adopt an average velocity of 41 km s~^ for the 
molecular material and obtain near and far kinematic distances of 2.7 and 10.5 kpc, using the 



Galactic rotation curve of Clemens (1985). Rathborne et al. (2006) chose the near distance for the 



IRDC based on the assumption that it must be closer than the bright IR emission in order to appear 
in silouhette. Also, the UKIDSS images exhibit a noticeable decline in source counts coincident 
with IRDCs relative to the field at large. Near-IR source catalogs from the UKIRT Infrared Deep 
Sky Survey (UKIDSS; Lawrence et al.||2007 ) show 80 - 160 sources per square arcminute coincident 



with the IRDCs, while the field away from the molecular shows typically shows over 400 sources per 
square arcminute. This indicates that a large fraction, up to 75% of sources are being obscured at 
the location of the molecular cloud. If the molecular cloud and associated IR features were at the 
far distance of 10.5 kpc, then the cloud would be behind a majority of the field stars and any decline 
in source counts would be much less noticeable. Therefore, we adopt the near distance of 2.7 kpc 
for the entire region. The difference between the lowest and highest peak emission velocities for the 
individual regions in Figure [2] amounts less than 0.2 kpc in distance, which is less than the assumed 



error in kinematic distances which we take to be 15%. The observed velocity spread is consistent 
with observed internal velocity dispersions of molecular clouds of several km s~^ (Solomon et al. 



1987 Ballesteros-Paredes et al.||2011 ). We estimate the column density of the G38.9-0.4 molecular 
cloud from the GRS ^^CO cube using the formula from Simon et al. (2001), 



N{^^CO) = 8.75 X 10^^ TmbAv cm" 



(1) 



where Tmb is the main beam brightness temperature in K and Av is the velocity range in km s ^, 
for which we used 35-45 km s^^. This equation assumes a ^^CO excitation temperature of 10 K, 



which is consistent with the estimated average value for the cloud of 8.3 K (Roman-Duval et al 

ICC 

45 and X{^'^C0/Il2) = 8 x 10"^ (|Simon et all200l[) and 



2010). We then compute the total gas mass surface density, Sco; in g cm ^ from ^^CO using 
the conversion factors R(^^CO/^^CO) 



multiplying by fJ-mn, the mean molecular weight (2.72; Roman-Duval et al. 2010) times the mass 
of neutral hydrogen. The total mass is the mass surface density summed over the physical size of 
an individual region, which was estimated using a distance of 2.7 kpc. The mass errors are from 
the flux estimation, the assumed distance, and the adopted temperature. Table [T] shows the results 
for the five main subregions and for the entire cloud. Column 1 gives the region name, columns 
2 and 3 list the center coordinates, and column 4 is the nominal diameter of each subregion. The 
specified areas are used to calculate the total gas mass for each subregion. Column 5 gives the peak 
velocity for each subregion and column 7 gives the ^^CO derived gas mass. 



2.2. Submillimeter Dust Continuum 

One of the key projects undertaken with the Herschel Space Observatory was the Herschel 



Infrared Galactic Plane Survey, or Hi-GAL (Molinari et al. 2010). Hi-GAL covers the Galactic 
Plane at 70, 160, 250, 350, and 500 /^m from | / |< 60° and | h |< 1°. The advantage of using 
Hi-GAL millimeter wavelength dust continuum over ^^CO line emission to trace molecular gas 
is that carbon-based molecules deplete onto dust grains when densities exceed 10^ cm~^, while 
mm-continuum remains optically thin at the higher densities expected in IRDCs and SF cores 
(Tafalla et al. 2004). While the ^^CO observations have the advantage of separating foreground 



and background emission to achieve a better estimate of total cloud mass, the sub- millimeter 
continuum should provide more accurate measurements in dense regions, such as IRDCs and SF 
clumps. 

The Hi-GAL wavelength coverage spans the peak of thermal dust emission (15 K at ~100 
;um), which allows for temperature fitting and more accurate surface density estimates. We made 
a dust temperature map by first aligning and resampling the Herschel maps to the 3'.'2 pixel scale 
of the 70 /xm maps using bi-linear interpolation. This approach has the advantage of perserving 
angular resolution and information in the short-wavelength maps but assumes that the emission at 
longer wavelengths varies smoothly over sub-pixel scales (14" pixel"^ at 500 /im). We then made 
integrated velocity maps of associated ^'^CO emission (35-45 km s~^) and the total line-of-sight 



emission up to the tangent point velocity (0-79 km s~^), and used the ratio of the two to determine 
the fraction of sub-mm emission associated with the SFR. This method makes the assumption that 
sub-milhmeter emission scales linearly with ^^CO emission, which may not be the case over the 
wide density range observed. The density where ^^CO becomes optically thick and where it depletes 
onto dust grains will alter the ratio. As a result of these two effects, we may be underestimating 
the associated sub-mm flux, leading to lower masses and surface densities. This is only important 
in higher density regions, such as IRDCs, and should not affect typical diffuse ISM. 

After subtracting possible contaminating foreground/background emission from the total in 



each band map, we used the IDL fitting routine MPFIT (Markwardt 2009) to fit the Herschel 



5-band data, pixel-by-pixel, to a modified blackbody of the form: 



lu = B^{T) {I - e"") MJysr 



-1 



(2) 



(3) 



where By{T) is the Planck function, and Tq is the optical depth at frequency Vq- For the fitting we 
chose Uo = (c/250 ^m), /3 = 1.75, and left t^, and T as free parameters. We allow the temperature, 
T, to range from 1 to 100 K and r25o from to 1. We also required data in at least three bands 



at any given pixel for be fit. Values for Vq and /3 were adopted from Rathborne et al. (2010) and 



Battersby et al. (2011). 



Figure |4] shows a grayscale image of the derived temperature map. Black contours outline 3 
K steps from 15-45 K, while green and magenta circles outline the different sub-regions as before. 
We derive average temperatures for off-cloud regions of 19-22 K, which agrees well with typical 



Galactic cirrus (Battersby et al. 2011). Fitted IRDC temperatures generally range from 14-22 K, 



and the bright rims of the bubbles have temperatures of 25-45 K. Computed average values for the 
subregions are listed in column 6 of Table [l| these temperatures fit with the general picture that 
IRDCs contain cold SF cores, while the bright-rimmed clouds are H II region shells being heated 
by the winds and UV photons from massive stars. 



Our procedures and results are similar to those from Rathborne et al. (2010), except they use 



a variable /?. We performed a trial fit with /3 variable from 0.5-2.5, and the final temperatures 
were generally within 20%. The largest changes were in the IRDCs where (5 dropped to 0.6-1.0 and 
temperature increased from 15 K to 20 K. The bright rims were also fit to lower /3 in the range of 
1.5, which increased the temperatures by 3-5 K. However, with at most 5 data photometric points, 
we deemed it necessary to reduce the number of free parameters and adopt a nominal fixed value 
for /3. 



From the temperature map, we calculated the pixel-by-pixel total gas mass surface density 



with 



^500 



lu 



By\l jKju^dust 



g cm 



(4) 



where l,y is the observed brightness at 500 /im in MJy sr ^, Ki, is the frequency-dependent opacity 
of 5.04 cm^ g~^, Xdust is the dust-to-gas ratio of 0.01, and B^(T) is the Planck function in MJy sr~^. 



The opacity was adopted from the tables of Ossenkopf & Henning ( 1994 ) assuming thin ice mantles 
at a density of 10® cm" 



■^ after 10^ yr. The surface densities are then multiplied by the area to 
arrive at a total gas mass, which for the IRDC is 1386ib418 Mq. Rathborne et al. (2006), using 



1.2 mm continuum data from the IRAM 30 meter and a dust temperature of 15 K, derived a mass 
of 1792 Mq, in rough agreement with our estimations from both ^^CO and 500 /im measurements. 
Full results are listed in column 8 of Table [Tj We calculated the surface density for each of the five 
bands and found similar results for 170, 250, 350, and 500 ^um, but 70 /xm varied significantly. We 
adopted the 500 fim values as this wavelength has the lowest optical depth and is likely to be the 
most reliable tracer of dense gas. 



2.3. Ionized gas as seen by radio continuum 

Figure [5] is a three-color image with 500 ;um in blue, 8 //m in green, and 24 /im in red. White 



contours outline the 1.4 GHz radio continuum from VLA Galactic Plane Survey (VGPS; Stil et 



al. 2006). The minimum level is the mean background temperature plus two times the standard 
deviation, which is 13.15 K + 2(0.43) K = 14.01 K. Subsequent contours are 4, 6, 8, 10 and 12a 
above the mean background temperature. The continuum images from the VGPS have 1' resolution. 
The figure shows strong sub-mm emission from cold dust coincident with the IRDCs and on the 
peripheries of the bubbles. The radio emission peaks inside the two bubbles N 74 and N 75 and 
exhibits a lesser peak near MCM 16. Fluxes were extracted for N 74, N 75, and MCM 16 using 
IMVIEW, but there was no detection for MCM 15. Considering, the weak emission from MCM 16, 
it may be appropriate to treat the flux as an upper limit. The radio emission within the bubbles is 
most likely free-free emission from H II, so we estimate the number of ionizing photons necessary 
to produce the 1.4 GHz flux by: 

-1 



Nl- 



•jc 



7.5 X 10^6 F, i.O-1 d^ r-°-45 



ph s 



(5) 



from Matsakis et al. (1976), where F^ is the flux density in Jy, u is the frequency in GHz, d is the 
distance in kpc, and Te is the electron temperature in 10^ K. We adopt 2.7 kpc for the distance. 



8800 K for the electron temperature ( Balser et al.||20lT >, and the frequency is of the VGPS is 1.42 
GHz. The results are summarized in Table [T| where column 9 is the measured flux, column 10 is 
the computed logio(NLyc) ph s^^ and column 11 is the 09.5V equivalent (NLyc(09.5V) = 10^^-^^ 
|2005| ). We find a UV flux of Niyc = 10^®°® ph s"^ for N 74, in agreement 
-|^q46.03 pj^ g-1 fpQj^ Beaumont & Williams (2010). Despite these 



ph s" 



-1. 



Martins et al. 



with the estimate of ^lyc 

values for the Lyman continuum flux, it is still likely that the bubbles N 74 and N 75 each contain 

a late-0 star (see Section |4]). Indeed, Watson et al. (2008) investigated the O star content of three 



prominent bubbles (NIO, N12, and N49 from Churchwell et al. (2006)) and found that the radio 
flux underestimated the ionizing flux by a factor of two when compared with estimates based on 
the known spectral type through optical spectroscopy. This is likely because a significant fraction 



of UV photons are being absorbed by dust and heating the grains, as evidenced by the bright 24 
lira emission within IR bright bubbles. Additional losses may be incurred if UV photons escape the 
bubbles without being absorbed; therefore radio estimations of Lyman continuum emission should 
be treated as lower limits. Using a single 09.5V star as the estimated O star content of N 74 and 
N 75, we can calculate the mechanical wind luminosity, L^, 



U 



0.5 rdw v:;^ erg s 



,2 



-1 



(6) 



where m^ is the mass-loss rate and v^ is the wind velocity. This yields IQ'^^-^* ergs s for nominal 
values of 10~^'^ M© yr~^ and 1500 km s~^ for the mass loss rate and wind velocity of an 09.5V 



star (Marcolino et al. 2009). We then calculate the H II region expansion radius as a function of 



time from Weaver et al. (1977) as. 



R 



27 L^ n-O-2 tO-6 



pc, 



(7) 



where L^ is the wind luminosity in units of 10^^ erg s~^, Uq is the ambient ISM number density in 
10^ cm~^, and tg is the expansion time in Myr. Figure|6jplots the radius versus age of expanding H II 
regions for three different initial average density values (10^, 10'^, and 10^ cm~^). The horizontal 
lines indicate the nominal radii of 0.65 pc and 1.1 pc for the N 75 and N 74 bubbles, respectively. 
The lower density (10^ cm~^) is typical of molecular clouds and higher values are approaching the 
density of IRDCs and cold molecular cores. For initial densities between 10^ ~ 10^ cm~^, we find a 
range of 0.4 - 0.8 Myr for the expansion age of N 75 and 1.0 - 2.2 Myr for N 74. RDI forms stars 
in a few 10^ yr after the ionization front has reached the existing clumps and CC may take a few 
Myr of H II region expansion for SF to begin. The apparent youth of N 75 may rule out SF via CC 
or RDI, while N 74 may be old enough for either process to have formed stars, depending upon the 
intial conditions 



3. YSO Identification and Properties 

We identified YSOs by fitting their observed near-IR (NIR) and mid-IR (MIR) spectral energy 



distributions (SEDs) to model SEDs. The procedure described below is similar to that of Robitaille 



et al. (2007) and Povich et al. (2009). 



3.1. Data Compilation 

We started with the point source catalogs (PSCs) from the Spitzer Space Telescope GLIMPSE 
I survey ( [Benjamin et al.||2003[ ). GLIMPSE I (hereafter, GLIMPSE) covered the Galactic plane 
(10° < I < 65°, I 6 |< 1°) with 1'.'2 pixel resolution in four bands, 3.6, 4.5, 5.8, and 8.0 /im. 
The GLIMPSE PSCs also include matched JHKg photometry from the Two Micron All Sky 



Survey (2MASS) (Skrutskie et al. 2006). We supplemented the 2MASS near-IR photometry with 



point source catalogs from the latest data release (DR7) of UKIDSS. UKIDSS has a finer angular 
resolution (C/S FWHM, 0'.'2 pixels at K) compared to 2MASS (1" FWHM, 2" pixels), as well as 
greater depth. The median 5o" limiting magnitudes for the UKIDSS Galactic plane survey (GPS) 
are J = 19.77, H = 19.00, and K = 18.05 ( [Lucas et al.][2008l ), while 2MASS is limited to J = 15.8, 
H = 15.1, and Kg = 14.3 (S/N = 10; Skrutskie et al. 2006). In addition, we performed aperture 



photometry on 24 ^m mosaic images (1'.'2 pixels) from the Spitzer MIPSGAL survey (Carey et 



al. 2009). The custom code calculated 24 /um fluxes and errors for the input source list using 



the location of GLIMPSE PSC detections and standard DAOPHOT-type routines. The relevant 
photometry parameters we adopted were an aperture of 7", an annulus of 7 - 13 ", an aperture 
correction of 2.05, and a zero-point flux of 7.17 Jy. If the extracted flux had a signal-to-noise ratio 
(S/N) less than five, then an upper limit, defined as five times the flux error, was set instead. 

UKIDSS matches to GLIMPSE sources were found by summing the astrometric errors (0'.'09 
and 0'.'3 for UKIDSS and GLIMPSE) in quadrature, which yields a matching radius of 0'.'313. If a 
UKIDSS match is found, then the UKIDSS photometry is used preferentially over 2MASS unless 
the magnitudes approach the saturation limit. The UKIDSS GPS has J = 13.25, H = 12.75, and 
K = 12.0 as saturation limits (Lucas et al. 2008), and we set limits 0.5 magnitudes fainter than 

13.75, H 



the quoted values. If a 2MASS magnitude is brighter than J 
then the 2MASS magnitude is used for that band. 



13.25, and K, = 12.5, 



The higher resolution and depth of UKIDSS images occasionally reveals single 2MASS/GLIMPSE 
objects to be multiple sources and increases the chance of a false match when trying to merge 
the data. If multiple UKIDSS objects are matched to a single GLIMPSE source, then the same 
GLIMPSE source is entered into the source list with different UKIDSS data points. However, for 
this study no multiply matched sources appear in the final YSO list. 

We require detections in a minimum of three bands to be included in the source list. Sources 
with three detections may still show a true IR excess even if their SEDs and properties are not 
well constrained. We fixed the minimum photometric error in any band to be 10% to prevent a 
single data point from dominating the resulting fits and to account for possible errors in absolute 



flux calibration (jRobitaille et al. 2007; Povich et al. 2009). Also, given the difficulty in extracting 



accurate MIPSGAL photometry, we use a minimum error of 15% for 24 /xm as in Povich et al. 



(2009). We did not make any color or magnitude cuts to our sources prior to SED fitting. 



3.2. Fitting Procedure 



Source SEDs from the initial list were fit to Kurucz (1993) stellar models using the routines 
of Robitaille et al. (2007). We allowed the extinction to vary from Ay = 0-35 magnitudes and 



picked x^/^data < 2 to be the threshold for a good stellar fit. Sources having SEDs well- fit by 
reddened stellar photosphere were removed. GLIMPSE sources with multiple UKIDSS matches 
were removed if one of them was well-fit to the stellar models. We also eliminated sources where 



10 



the IR excess appeared to be caused by a single band with large photometric errors. 

Sources that could not be matched to reddened stellar photosphere were then fit against a grid 



of 200,000 YSO model SEDs (Robitaille et al. 2006) that were computed using radiation transfer 



code presented in Whitney et al. (2003). YSO fitting used the same toolkit as stellar fitting 
( [Robitaille et al.||2007[ ). We chose model aperture sizes of 3" for 2MASS and GLIMPSE bands, 
l'.'5 for UKIDSS, and 7" for 24 //m and set the distance range to 2.4 to 3.0 kpc. The aperture and 
distance range eliminates models with envelopes that would be extended and/or partially resolved 
at the estimated distance to the region. The interstellar extinction was allowed to vary from Ay 
= 5-40 magnitudes, where the range was estimated from the minimum and maximum extinction 
values estimated from the Herschel maps. After fitting, we verified the IR excess by-eye in the 
images and SED plots. In total, we find 162 YSOs candidates in the G38.9-0.4 field. Table [2] 
lists the YSO ID number we assigned, the GLIMPSE catalog ID, Galactic longitude and latitude, 
and the observed magnitudes. The J, H, and K columns represent both 2MASS and UKIDSS 



magntidues according the selection rules from Section 3.1 



3.3. Physical YSO Properties 

The YSO models are calculated over a range of properties including, but not limited to: the 
total bolometric luminosity, orientation, mass of the central star, disk mass, and envelope accretion 
rate. Since each SED can generally be fit to many YSO models equally well, it is necessary 
to calculate the average YSO properties for a range of well-fit models. We chose models where 
(x^ ~ Xbest)/''^data < 1 to be acceptable fits. The weighted average was calculated by: 

P, = e-^^^#^, (8) 

UX^-A,yP, N 

^^ = E^F^ ^ W^V ^^^) 

where Pj is the probability relative to the best fit model, X is the quantity being averaged. Ax 
is the weighted average, N is the number of acceptable model fits, and o"^ is the variance in the 
weighted average. Table [3] gives values a few important YSO properties. Column 1 is the YSO 
index; column 2 is an number where each digit represents a magnitude (column) from Table l2| A 
'0' indicates that the band was not used in the fit, a '1' indicates the band was used, and '3' means 
the band was used as an upper limit, not a true detection. Column 3 is the number of acceptable 
fits, and column 4 is the x^ of the best fit model. Columns 5 through 12 give the best fit and 
average values of the stellar mass, disk mass, envelope accretion rate, and integer stage. When we 
refer to YSO stages for individual sources, we use the stage given in column 12. Columns 13 and 14 
are the average fractional stage (F stage) and dispersion in the acceptable fitted stages {a stage)- The 
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stage classification is based off the criteria givein by Robitaille et al. (2006). They assign stage O/I 
(hereafter, I) for heavily embedded sources with an envelope accretion rate M jM^ > 10~^ yr~^. 
Stage II YSOs have MjM^. < W'^ yr"! and Mdisk/M^ > 10"^. These objects have accreted 
most of their mass, but are young enough to still have a thick disk. Thin disk, stage III, YSOs have 
only a moderate IR excess and are characterized by M/M^ < 10~® yr~^ and Mdisk/^* < 10~^. 
Since YSO evolution is a continuous process, we have calculated the average fractional stage (Fstage) 
and uncertainty (astage) in addition to the integer stage. Fstage and cTstage give a representation of 
which models fit the object and how well constrained they are. In some cases agtage = 0; this 
does not imply a perfect fit. Rather it means that all of the well- fit models for that object give the 
same stage. This typically occurs when there are a limited number of acceptable fits. We calculate 
the integer stage by requiring that at least two thirds of the total weight (probability) be matched 
with a single stage: 

Emotai) > ^-^^ ' ^''^ 

so, if the combined weight of all stage I fits for a given object is greater than 0.67, that object 
is classified as stage I. If no stage has a combined weight greater than 0.67, then that object is 
classified as ambiguous. Using this method, we find 54 stage I, 70 stage II, 6 stage III, and 32 
ambiguous sources. Figure [7] shows an example SED for each of the YSO stages. The upper-left 
plot shows the SED for Y133 (stage I), the upper-right shows Y112 (stage II), the lower- left shows 
Y157 (stage III), and the lower-right shows an Y30 (ambiguous). The solid black curve is the full 
SED best-fit YSO model, the grey curves other well-fit model SEDs with (x^ — Xbest)/^data < 1, 
and the long-dashed curve is the SED for the underlying star if there were no IR excess. There are 
a wide range of well-fit SEDs for the ambiguous sources, which explains the difficulty in obtaining 
a stage I or stage ii classification. The plots clearly show the importance of far-IR photometry in 
constraining YSO stages. 

It is important to note that while the model SED fitter can give estimates of various YSO 
physical parameters, the estimates become more uncertain as the number of data points and wave- 



length coverage decreases. Deharveng et al. (2012) showed that the SED fitter may significantly 



overestimate the mass of heavily embedded YSOs with large envelopes, such as the stage I YSOs 
present in the IRDCs. It is also possible that there undetected multiple sources, so luminosities 
and mass may be overestimated in these cases. Thus, the physical parameters given in Table [3] 
should be considered carefully. 



3.4. G38. 9-0.4 YSOs 

Color-color and color-magnitude diagrams (CCDs and CMDs) are also useful for identifying 
and classifying YSOs, and there are several different schemes employed (Allen et al.|2004 Rho et al 



2006t [Gutermuth et al.|r2008D . Figure [s] plots mid-IR CCDs and CMDs of YSOs identified through 



their full SEDs. YSOs are shown as plusses for stage I, triangles for stage II, squares for stage III, 
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and circles for ambiguous YSOs. Grey symbols indicate that the only have an upper limit at 24 
^m. The arrows show an extinction of Ay = 15 magnitudes (Indebetouw et al. 2005) The colors 



and magnitudes were chosen to provide the greatest separation between YSOs and main-sequence 
(MS) stars and between different YSO stages. The [5.8] — [8.0] color shows excellent separation 
between YSOs and the MS locus in both the CMD and CCD. The YSO stages do not separate well, 
although the stage I objects are more red in [3.6] — [4.5] than stage II or ambiguous YSOs. The 
[3.6] — [4.5] versus [8.0] — [24] CCD is much better at separating stage I YSOs from stage II YSOs. 
It is interesting to note that in each of the four diagrams the ambiguous objects tend to fill the 
gaps in between stage I and stage II YSOs. This may indicate that they are in a transition phase 
between classical stage I and stage II. This is also evident in the average properties in Table |3j 
Four YSOs stand out with [5.8] - [8.0] < 0.25 (Y36, YSO, Y117, and Y130). Y36 has a weak [5.8] 
detection and Y117 has a weak 8 /im detection, which is the likely cause of the blue color. Y50 
shows a broad dip in flux near 9 /xm, which may be a result of self-absorption from viewing an 



edge-on disk (Whitney et al. 2003). The fits to Y130 show a more narrow dip in flux, which may 



be caused by silicates at 9.8 /im (Whitney et al. 2003). 



There are several possible contaminants in any YSO sample, the most common of which are 



asymptotic giant branch (AGB) stars (Robitaille et al. 2008). AGB stars represent an advanced 



stage in the evolution of low- and intermediate- mass stars (Herwig|2005>. They are also intrinsically 



cool and are prone to mass loss resulting in significant dust formation (Groenewegen 2006), and 



so they radiate strongly in the IR with YSO-like near- and mid-IR colors (Whitney et al. 2008). 



Robitaille et al. (2008) published a catalog of intrinsically red sources for GLIMPSE I/II, which 



includes preliminary classifications for YSOs, AGB stars, and unresolved planetary nebulae and 
galaxies. One of the YSOs in the initial list was listed as an AGB star in the catalog and was 



removed from the final YSO list. Whitney et al. (2008) found that most AGB stars could be 
removed by eliminating sources with [8.0] - [24]<2.2. There are 131 YSOs with both 8 and 24 fim 
detections and only three (2%) fall below this color criterion. In addition, AGB stars are also likely 



to outnumber YSOs when considering sources with [4.5] < 8 (Povich et al. 2009); only one YSO is 
brighter than this threshold as the brightest sources typically presented SEDs that were removed 
earlier in the YSO selection process. 

In order to estimate possible contamination of our list by unassociated field YSOs we performed 
the same analysis on a nearby, quiescent field region. It was a circular region equal in area to the 
target region with a radius of 12'.'3 centered on / = 38? 4490 and b = -0?4548. For G38. 9-0.4 we 
found 339 sources poorly fit to reddened stellar photospheres with 162 being bona fide YSOs. In the 
field region, we find 146 poorly fit sources with just eight candidate YSOs. The G38.9-0.4 region 
contains 176 contaminating sources (not MS stars or YSOs) and the field region contains a similar 
number, 138. The 28% difference is most likely due to natural variation in the field star density, 
but a fraction of the extra sources might be YSOs that were removed due to weak IR excess or 
unreliable photometry. The ratio of field YSOs to SFR YSOs is quite small, 5%, so we can assume 
that most of the YSOs in Tables [2] and [3] are associated with the SFR itself. We also find that there 
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are no AGB stars in the field as identified in the Robitaille et al. (2008) catalog while there was 



just one in the SFR, so we also conclude that the AGB contamination towards G38.9-0.4 is small. 



Povich et al. (2009) estimate < 20% contamination in their YSO list for M17, however, our region 
is covers a much smaller area, is further from the Galactic center (l ~ 15° for M17), and likely has 
a less complicated line of sight. Therefore, we may expect to have a lower contamination fraction. 

Figure [9] is a three-color image of G38.9-0.4 with UKIDSS K, 8 //m, and 24 fim in blue, green, 
and red, respectively. Red pluses, green triangles and cyan circles indicate stage I, stage II, and 
ambiguous YSOs. White contours delineate the gas surface density levels of 0.01(1.75*), equaling 
0.01, 0.0175, 0.0306, 0.0536, 0.0938, and 0.1641 g cm'^. The YSOs tend to cluster near the IRDCs, 
with several stage I sources surrounding MCM 15. Most of the YSOs around N 74 are stage II or 
ambiguous, which indicates that it is older than MCM 15 and, likely, older than N 75 as well. The 
high level of diffuse 8 and 24 //m emission coincident with MCM 15, which may be obscuring YSOs 
and causing a higher level of incompleteness. Given the apparent young age and signs of recent SF 
in this region, we might expect to find more YSOs coincident with the bright bubble rims, but this 
is not seen. It is possible that YSOs have not yet, or will not be, formed in the bubbles shells, or 
they are simply hidden amid the high IR background. 



3.5. Comparison to Other Regions 

The final YSO count, 54 stage I, 70 stage II, 6 stage III, and 32 ambiguous, is quite different 
than those in other known SFRs. The observed ratio of stage I YSOs to stage II of 0.76 is quite 
high compared with other reported values of 0.55, 0.4, 0.38, and 0.01 ( Kang et al.|2009 Gutermuth 



et al. 2008 Povich et al. 2009, Koenig et al. 2008 respectively). Each of the regions studied are 



very different in terms of morphology, sky coverage, and YSO identification methods, which may 
be responsible for some of the differences. It is also possible that this represents a possible age 



difference with higher stage I/II ratios representing younger regions. Povich et al. ( 2009 ) studied 
M17 and Koenig et al. (2008) studied the W5 H II region, both of which have distances of ~ 2 kpc. 



while the main clusters have ages of 3 Myr (M17) and 5 Myr (W5) based on the maximum ages 
of the most massive stars. These ages, and the 1-2 Myr ages for N 74 and 75, are significantly 



longer than the average lifetime of stage I YSOs (0.1-0.2 Myr; Andre et al. 2000). Both M17 and 
W5 show evidence for multiple epochs of SF ( [Povich et al. 2009 Koenig et al. 2008), where the 
initial populations may have evolved beyond the embedded stage I phase lowering the overall stage 
I/II ratio. However, the G38.9-0.4 SFR may still be in the first generation of SF and thus have a 
younger YSO population. The large fraction of YSOs coincident with IRDCs, which are sites of 
current SF, is also suggestive of young stars. The lack of stage III objects, and stage II to some 
extent, may be explained by the difficulties in detecting weak IR excesses, especially in regions of 
high and variable extinction. A stage III YSO SED may show little-to-no excess and appear stellar 
shortward of 10 /im, and these may be removed by the initial stellar SED fitting. Stage II YSOs 
can show excesses at shorter wavelengths, but extinction may hinder their detection. 



- 14- 

4. Breaking Down G38. 9-0.4 

4.1. Field Star Identification 

The location of G38.9-0.4 in the heart of the Galactic Plane makes it difficult to assess the 
associated stellar content, in part because of extinction that is high and highly variable. Also, field 
stars likely outnumber associated MS and pre-main sequence (PMS) stars and must be removed in 
order to recover likely cluster stars. In order to look for stellar clusters we follow the prescription 



of Maia et al. (2010) to evaluate and remove the field star contamination. 



The field star estimation tool compares a target region to a field region using a UKIDSS 
JHK 3-dimensional color-color-magnitude diagrams (CCMD). After investigating different color 



and magnitude combinations, we used the J, J — H, and J — K magnitudes, as in Maia et al 



(2010). CCMDs were constructed for both the target regions containing a putative cluster (about 
10" radius for MCM 15, MCM 16, N 75, and N 74) and a field annulus around the target center. 
We made the annulus as large as possible without overlapping a neighboring bubble, 1-2' for the 
outer radius. Since the target and field areas are, in general, not equal, the number of field stars is 
normalized by the target-to-field area ratio, A. 

To estimate the field star content, the CCMD is divided to 3D bins along the magnitude and 
color directions. The cluster membership probability is then determined by the ratio of stars in 
the target bin stars to those in the field bin, P = (Nt — Nf x A)/Nt. So, if there are five stars in 
the target bin and one star in the field bin, then one star in the target bin is a "field star" and 
the membership probability for each target star is 0.8. In cases where the number of field stars 
in a cell is equal or greater than the number of target stars in the same CCMD cell, membership 
probabilities are zero. 

In order to mitigate sensitivity to the choice of bin center and bin size, we used three different 
cell sizes for each color and magnitude. Bin widths were 0.15, 0.30, and 0.45 mag for J — H and 
J — K, and 0.5, 0.75, and 1.0 mag for J. For each different bin size, we shift the cell centers in 
each dimension by one-third of the cell width in each direction. There are a total of six parameters 
(three width dimensions and three shift dimensions) each with three different values, this yields 
36 = 729 different grid permutations. For each permutation, the stars in the target cell are assigned 
a membership probability described above. They are averaged over all iterations to arrive at a final 
membership probability for each star in the target region. We tested our code on cluster LKIO and 
found that we could recover the morphology of the cluster CMD shown in Figure 5 of Bonatto &; 



Bica (2009). 



4.2. MCM 15 



Figure 10 shows a three-color UKIDSS JHK image of MCM 15 (left panel). White contours 



outline MIPS 24 /xm emission, and the magenta circle denotes the aperture used as the target region 
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for finding clusters. YSOs are shown as red plusses for stage I, green triangles for stage II and cyan 
circles for ambiguous sources. In the image we see YSOs on the outskirts of the cluster, situated 
over regions of high gas surface density (Figure p|. There are only one or two YSOs detected in the 
central area of the cluster, however, the IR background is extremely high, as evidenced by the 24 
fiiQ contours, so the data likely suffer from incompleteness. There appears to be an over density of 
stars, suggestive of a stellar cluster, accompanied by diffuse K-band emission, probably reflection 
nebulosity or emission from hot dust. There is a strong morphological similarity between the diffuse 
K-band and 8 /im emission where PAH emission is the dominant source of luminosity (Figure ^ . 



Figure 10 (right panel) shows J versus J — K CMDs (top) and the J — H versus J — K CCDs 
(bottom) for the full target area (left), the field region (center), and the cleaned target (right). 
The cleaned CMD and CCD show only sources with a membership probability > 2/3. The filled 
diamond and asterisk mark the locations of an AO and an 09.5V star along a ZAMS ischrone (solid 
line; Marigo et al. 2008). The dotted line curve marks a 1 Myr PMS isochrones (Siess et al.||2000). 



All isochrones have been reddened by Ay = 12, according to Cardelli et al. (1989), and placed at 
the assumed distance to the molecular cloud (2.7 kpc). The solid arrow shows the magnitude and 
direction of Ay = 5 magnitudes. The CCD shows an over density of stars having high membership 
probability near J — K = 2.5-3.0, which implies approximately 12 magnitudes of extinction for a 
population of MS stars. The CMD shows an excess of stars between the ZAMS and PMS isochrones, 
which we interpret as additional evidence of a young stellar cluster. The candidate cluster members 
show a significant dispersion in color, consistent with a high degree of differential reddening. Some 
of the brightest cluster stars with high membership probability lie near the PMS isochrones, and 
these roughly correspond to the luminosity of mid-B stars. This is consistent with a recent study 
suggesting that MCM 15, which is coincident with IRAS 19012+0505, is an intermediate-mass star 



forming region that will not form stars earlier than a spectral type B3 (Arvidsson et al. 2010). 



4.3. N 75 



Figure 11 is a three-color image of N 75 (left panel) with color diagrams (right panel) as 
before. Colors and symbols are the same as Figure [Toj There is a possible compact cluster within 
the 15" (0.2 pc) radius magenta circle, centered on the bright star 2MASS 19034727+0509409 (J 
= 11.35, H = 10.61, K = 10.13) marked by a yellow 'x'. The contours showing 24 /im emission 
exhibit a local maximum at this location, while the Herschel dust column density maps show a 
local minimum inside the bubble (Figure [S]). This localized peak seen at 24 /im is not visible in in 
the Herschel maps at 70 /im, leading us to interpret this feature as emission from very small dust 
grains heated locally by the star (or stars) within the compact cluster. Interestingly, the 1.4 GHz 
radio continuum (Figure pi) also peaks at this location. 

The CMD (top right panel) shows a sequence of ~12-15 stars that is well-matched to the PMS 
isochrones (2.7 kpc distance; Ay = 8) and a small grouping of ~7 extremely red stars, possibly 
background giants. This may indicate the presence of a young stellar cluster where the massive 
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OB star has arrived on the MS and others are still in the disk-bearing PMS phase. The color-color 
diagram (lower right panel) shows a grouping of stars near J — K = 2.0, consistent with a stellar 
cluster having similar reddening, in this case A\i = 8 mag. The bright star, J = 11, J — K = 1.1 
on the CMD, is consistent with a 09.5V star at 2.7 kpc with ^v = 8 mag. Figure [12] plots the 
Teff — R relation for Kurucz stellar models that are well-fit to the observed SED (black dots) at 
the adopted distance of 2.7 kpc. The observed Teff — R relations for O stars from [Martins et al 



(2005) are plotted as solid gray lines and the location of 09.5V stars are shown as asterisks for 



the different luminosity classes. Where these two relations cross provides a good indication of the 



stellar spectral type (Watson et al. 2008), in this case a late-0 type star. Its location near the 



symmetry axis of the MIR nebula (partially seen in Figure 1 1 ) also suggests that this star may be 



responsible for creating the prominent ~1 pc diameter IR bubble (Figure [T]). Although the angular 
resolution of the VGPS is 1', the coincidence of the radio continuum peak with the 24 /um peak 
and the candidate 09.5V star implicates this star as the source of ionizing photons responsible 
for generating the (presumably) free-free emission filling the bubble. The remainder of the stars 
in the CMD along the PMS isochrone are consistent with mid-B to early-A stars. As such, their 
mechanical luminosity and ionizing flux would be negligible compared to a 09.5V star. Hence, 
the N 75 bubble appears to surround a compact young cluster dominated by a single late-0 star, 
but consisting primarily of intermediate-mass (and unseen lower-mass) stars. This cluster may be 



a descendant of an ultracompact embedded cluster (Alexander &; Kobulnicky 2012) and a more 



massive version of a Herbig Ae/Be cluster (Testi et al. 1997, 1999). 



4.4. N 74 



Figure 13 shows the three-color JHK image (left panel) and color diagrams (right panel) for 
N 74; colors and symbols are the same as Figure [TO} The white 24 /im contours peak over another 
potential compact cluster, anchored by the bright star 2MASS 19035684+0506370 in striking sim- 
ilarity to the N 75 bubble. Figure [5] also shows a local maximum in radio continuum flux near this 
cluster. 

The CMD (top right panel) shows a single bright star (J = 11.21, H = 10.35, K = 10.01) 
plus a sequence of stars that is well-matched to the PMS isochrone (2.7 kpc; Ky =8). This is 
again suggestive of a stellar cluster with properties similar to that of N 75, a single 09.5V that has 



evolved onto the MS while its lower mass companions are still in the PMS phase. Figure 14 shows 
the Teff — R relation for model fits as in Figure 12 , but for the candidate ionizing star in N 74. The 
data are consistent with this star having a late-0 spectral type. 
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4.5. MCM 16 



Figure 15 shows the three-color JHK image (left panel) and color diagrams (right panel) for 



MCM 16; colors and symbols are the same as Figure 10 The UKIDSS JHK image (Figure [15| left 



panel) again shows peaks of 24 /im emission and radio continuum (Figure [5]) centered near a bright 
star, 2MASS19041919+040641 (J = 11.43, H = 10.78, K = 10.48), however, the cluster is not as 
compact as in N 75 and N 74. This region is further away from the main molecular complex, so 
there is far more stellar background contamination, as evidenced by the apparent increase in source 
density. This region also lacks the dense molecular gas and YSOs that are present in MCM 15, 
N 75, and N 74. This region is likely to be the oldest in the complex based on the extent of the 
mid-IR emission (Figure [I|. The shape of the 8 and 24 ^m emission does not provide a clear radius, 
however, we estimate the age to be 2-3 Myr based on the closest IR rim to the cluster. This more 
advanced age agrees with the cluster being less compact as it has passed through the gas dispersal 
phase and may be dissolving. 



The CMD (Figure 15, top-right panel) shows a similar structure to the previous clusters, 
but with a lower extinction (^v = 6.5) and some additional scatter, likely due to the increased 
contamination by field stars. There is another bright star with colors consistent with an 09.5V, 
however, the Tgff — R relation shows this star to be inconsistent with a late-0, and it may instead 
be an early-B type star. 



5. Discussion 



5.1. YSO Mass Function 



Figure 16 plots the cumulative YSO mass function (YMF). The solid black histogram is the 
observed YMF, the solid grey line is the best fit slope of Vymf = —3.8 it 0.1, and the dashed 

We have chosen to use a cumulative 



grey line is the Salpeter (1955) slope of r555 



-1.35. 



histogram to mitigate the effects of bin size and bin location, which become important with small 
sample sizes (i.e., high mass bins). In addition, the slope estimation was restricted to masses > 2.5 
Mq to account for lack of completeness at lower masses. The resulting best fit slope is much 



steeper than the Salpeter IMF, but agrees quite well with the YMF found by Povich et al. (2011) 
of Vymf = -3.2 for the Carina Nebula and Vymf = -3.5 for M17 SWex (iPovich & Whitney 



2010). The YMF slope for our sample is likely dominated by incompleteness, especially below 



~ 2 Mq where source counts drop off, but inclusion of such sources would only make the fitted 
slope steeper. On the more massive side, the lack of higher mass YSOs (> 6 Mq) may be due to 
stochastic sampling of the IMF or possibly that high-mass stars evolve to the MS more rapidly than 



low-mass stars (Povich &: Whitney 2010). Figures 11, 13, and 15 each show a likely star cluster 



with a massive central star. The most massive star appears to be on, or close to, the MS while the 
fainter stars seem to fall along the PMS tracks. 
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The most likely scenario for explaining the lack of massive YSOs is simply that this SFR does 



not present the appropriate conditions. Krumholz & McKee (2008) find a minimum surface density 
of 1 g cm~^ for the formation of massive stars, while the highest surface density in the Herschel 500 
/xm map is ~ 0.25 g cm~^, and values above 0.1 g cm~^ are generally limited to the IRDCs. This 
suggests that, at present, this SFR does not have densities high enough to form massive stars. This 
suggestion is supported by the classification of the MCM 15/N 75 region as an intermediate-mass 



SFR (Arvidsson et al. 2010). 



5.2. YSOs and Gas Surface Density 



Figure 17 shows a plot of individual YSO mass versus gas mass surface density derived from 
Herschel 500 //m maps. The upper left plot shows stage I YSOs, upper right shows stage II, lower 
left shows ambiguous, and lower right shows stage III. The heavy blue symbols are those sources 
within the radii of the N 75 and N 74 bubbles given in Table [T| and extra large symbols show the 
arithmetic average of all the sources in each panel. Since it takes more gas to form a larger star, one 
might expect more massive YSOs to be embedded in higher surface density regions than less massive 
YSOs, however, this trend is not seen in the plots. The linear correlation coefficients are 0.22, 0.19, 
and 0.22 for stage I, stage II, and ambiguous YSOs, which indicates a weak degree of correlation 
for 54, 70, and 32 data points, respectively. There might be evidence for an upper envelope as 4 
Mq YSOs are typically associated with a surface density higher than 50 Mq pc~^ (0.01 g cm~^). 
The estimated uncertainty on an individual YSO mass ranges from 0.2-0.5 Mq depending on how 
well constrained the models are, but variations of this size appear insufficient to obscure any tight 
underlying correlation. However, there are likely to be other sources of scatter in the estimated 
masses and gas column densities. The difference in angular resolution between the GLIMPSE and 
Herschel 500 /im beamsize (1'.'2 vs. 36" FWHM) contributes uncertainty in matching a given YSO 
with a local gas mass surface density. The surface density maps were interpolated to the 3'.' 2 scale 
of the 70 ^m images, but this is not equivalent to sub-mm data with a small native pixel scale. 
If the ratio of unrelated background gas to G38.9-0.4 molecular cloud gas varies on scales smaller 
than the GRS beam size, then additional pixel-to-pixel noise will be added to the surface density 
maps. However, we do find that, on average, stage I YSOs fall at a higher surface density (~239 
Mq pc~^) not quite twice that of stage II YSOs (^141 Mq pc~^). The ambiguous sources fall 
in between these values (~191 Mq pc~^), which may indicate that these objects are transitioning 
from a heavily embedded phase to a less embedded phase. While there are few stage III YSOs, 
it is worth noting that they fall in the lowest density regions with an average gas density of ~ 112 

Mq pC-2. 



Figure 18 is a plot of YSO mass surface density against the local ISM gas mass surface density. 



We followed Gutermuth et al. (2009) and calculated the angular radius out to the sixth nearest 
neighbor YSO. Then, using the adopted distance of 2.7 kpc, we converted the angular radius into 
a physical radius and calculated the enclosed area. The YSO mass surface density is then the 



19 



sum of the masses of the N=6 YSOs divided by the physical area they cover. The upper left 
panel shows all of the YSOs, upper right shows stage I YSOs, lower left shows stage II, and lower 
right shows ambiguous sources. Heavy blue symbols are those YSOs within the N 74 and N 75 
boundaries. We find excellent power-law correlations between the YSO mass densities and gas 
surface densities with the best fit equation given in the upper left corner of each panel. The slopes 
are 1.27±0.06, 1.38±0.11, 1.33±0.08, and 1.30±0.17 for ah YSOs, stage I, stage II, and ambiguous 



sources, respectively, so within the uncertainties, the slopes are equivalent. Gutermuth et al. 



(2011 ) find a wide range of slopes between 1.4 and 3.8 for eight nearby SF regions, but they did not 



calculate slopes for distinct evolutionary stages. They also adopt an average mass for their YSOs 
of 0.5 Mq, while we use the average individual mass from the model YSO fits. We also use the 



sixth nearest neighbors, while Gutermuth et al. (2011) use up to the 11th. If the nearest neighbor 



count was increased from 6 to 11, then the resulting slopes all drop by about 11%, meaning that 
our results are relatively insensitive to the choice of N. Another difference is that the regions in 



studied in Gutermuth et al. (2011) have distances less than 1 kpc and thus have a greater ability 



to detect low-mass YSOs. Figure [T6| shows a clear drop in YSO counts below 2 Mq, which may 
have significant effects on the slope of the YSO-gas density relation. If undetected YSOs have a 
similar spatial distribution to the current YSO list, then the slope of the density relation would 
remain the same but with a larger y-intercept. However, it is possible that faint (i.e., low- mass) 
YSOs are preferentially detected in areas with lower extinction and lower sky background. If this 
is true, then a higher fraction of low-mass YSOs may remain undetected in high gas surface density 



regions. This type of detection bias could steepen the slope shown in Figure 18 and bring the values 



more in line with those from Gutermuth et al. (2011). Interestingly, the Schmidt-Kennicutt Law 



(Kennicutt 1998), which describes the power-law relation between the areal SF rate density and 



the mass surface density has a power law index of 1.4 similar to the YSOs in this region. 



5.3. Evidence for Triggered SF? 

In the Introduction, we described the basic mechanism for three theoretically defined triggered 
SF scenarios: CC, RDI, and RR. Another mode of triggered SF involves the sequential collapse 



of material along a molecular filament. Simulations by Fukuda &: Hanawa (2000) show that as an 



H II region expands in a filament, it will cause the cloud to "pinch" and trigger the collapse of 
new SF cores to either side of the original H II region. If these new cores eventually form their 
own H II regions, the process will repeat. Recent evidence has shown that SF molecular clouds 



tend to be filamentary in nature (Molinari et al. 2010), and observational studies have suggested 



that "filament pinching" might trigger sequential waves of SF (Samal et al. 2012) and work by 



Koenig et al. (2012) hints at a similar sequential SF scenario. In regions with clear MS populations 



and known massive stars, they argue that the distribution of YSOs mirrors the pre-SF filamentary 
structure of the parent molecular clouds. They also suggest that massive stars and H II regions 
enhance SF by compressing the filaments. 
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The GLIMPSE 8 fira and Herschel 500 nva. images (Figure ^ show that the spaces between 
the bubbles have dark lanes. This may indicate a swept-up shell, or, it may instead be evidence 
of a pre-existing filament that extended from north of MCM 15 through the bubbles and further 
south. In this model one of the bubbles could have been the initial H II region, which triggered 
the sequential collapse of the cloud. However, the estimated bubble ages are similar and it is not 
clear whether one would have had sufficient time to trigger the other. It may be that the bubbles 
formed within the same cloud but do not have a causal relationship. 



In a large statistical study, Thompson et al. (2012) cross-matched the location of massive YSO 



(MYSO) candidates ( [Urquhart et al.|2008D with IR bubbles ( [Churchwell et al.|2006D and found that 



22% of MYSOs are coincident with bubble rims. Kendrew et al. (2012) used the more complete 



bubble catalog from the Milky Way Project (Simpson et al.|[2012) and found a similar coincidence 



within the bubbles. They also found that 67% of MYSOs lie within 2 bubble radii. These data 
suggest that triggered SF occurs in the presence of bubbles and H II regions. However, many 
bubbles are likely to reside in or near large SFRs and molecular cloud complexes. As such, bubbles 
are not unbiased regions distributed at random, but rather they are likely biased toward higher- 
than-average column density regions where stars are already expected to be forming. Figure [9] 
shows that most of the YSOs may be associated with the IRDCs surrounding N 74 and N 75 rather 
than with the bubbles, even at distances exceeding one bubble radius. Most of the YSOs that lie 
between one and two bubble radii away from N 74 and N 75 are most likely associated with nearby 
IRDCs or other features and the bubbles. Therefore, extending the triggered region beyond one 
bubble radius increases the likelihood of including spurious YSOs that are not causally connected 
to the bubbles. 

One method of demonstrating the occurrence of triggered SF is to establish a radial YSO stage 
gradient as a proxy for age. The N 74 bubble has two stage I YSOs within it compared to six stage 
II YSOs, while the rim is a mixture of stage I, II, and ambiguous YSOs. The bubble interior then 
has a stage I/II ratio of 0.2 compared to the 0.75 stage I/II ratio seen over the entirety of G38.9-0.4. 
This may be suggestive of an age gradient, but the limited number of YSOs around the N 74 and 
N 75 prohibit a statistically significant determination. We also see a dip in gas surface density 
(Figure [5]) coincident with the bubble centers and a density enhancement on the periphery, which 
is suggestive of a swept up shell. However, the portion of the N 74 rim that shows the largest 
number of YSOs also abuts the IRDC, and the YSOs may have formed from the IRDC prior to 
the advancement of the H II region. We conclude that there is no definitive evidence for an age 
gradient in the N 74/N 75 bubbles. 

Triggered SF via RDI and RR is characterized by a relatively short timescale (< 1 Myr) and the 
presence of YSOs with trailing, wind-blown pillars pointed towards the exciting source(s). There 
appear to be no such structures within the entire region. It is possible that the flux of ionizing 
Lyman continuum radiation is not high enough to effectively compress the surrounding ISM to 
form stars or that not enough time has passed for YSOs to have formed through RDI. However, 
at the estimated bubble age of 1 Myr, RDI triggered YSOs should have had enough time to form. 
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We conclude there is no evidence that RDI plays an important role in this SFR. 

Triggering via CC predicts that SF will occur in massive, regularly spaced clumps over Myr 
timescales and that the subsequent generation will include massive stars (Elmegreen &: Lada||1977). 



N 74 and N 75 appear to be wind-blown bubbles and the Herschel sub-mm and gas surface density 
maps indicate that they have been at least partially evacuated, as would be seen in a CC triggered 
SF scenario. However, there is no evidence for current SF within the bubble rims, and only a 
handful of YSOs have been detected within the bubbles, none of them massive. Also, while the 
gas surface density does drop within the bubbles it is not obvious that gas condensations have 
formed as a result of stellar feedback. The highest gas surface densities detected near the bubble 
peripheries (0.05 g cm~^at the upper left corner of N 75 as seen in Figure [I]) occur in the same 
azimuthal direction as the IRDCs that lie just outside the bubbles. The projection of the IRDCs 
and bubble rims makes it difficult to determine whether the gas is swept-up material from the 
bubble interior or the ambient molecular gas comprising the IRDCs, or both. We suspect that the 
majority of the gas is related to the IRDCs because surface density contours match the outline of 
the IRDCs visible in the GLIMPSE 8 ^va. maps. For N 74 and N 75, we conclude that there is no 
compelling evidence to invoke triggered SF through CC on bubble rims as a result of the stellar 
feedback. 



5.4. Does stellar feedback enhance or suppress SF? 

If feedback enhances SF through triggering, one possible observable consequence is a higher 
star formation efficiency (SFE), i.e., a higher density of YSOs per unit gas mass in triggered SFRs 
than in non-triggered regions. That is, triggering may elevate the SFE by precipitating collapse of 
clouds that are marginally stable and otherwise would not have formed stars. Numerical simulations 
suggest that stellar feedback may both trigger and suppress SF simultaneously and not have a strong 



effect on the overall SFE (Dale et al. 2012). Our analysis is consistent with a scenario whereby 
stellar feedback does not have an appreciable net effect on the total SFE of a molecular cloud. 
Figure [18] does not show any significant difference in the density of YSOs on the rims or within 
bubbles (blue symbols) compared to those not obviously influenced by massive star feedback (black 
symbols; regions outside the bubbles). Although the statistical power in this sample is limited, 
this null result may he taken to mean that triggering either has not occurred, or that feedback does 
not affect the SFE and is not as significant a factor as the initial molecular cloud density. This 
conclusion is in agreement with others who have suggested that triggering does not measurably 



affect the global SF rate (Elmegreen 2011). Stellar feedback-induced triggering may be just one 
of many mechanisms to form stars among processes such as the isolated spontaneous collapse of 
molecular cores, SF driven by cloud-cloud collisions, or other dynamical effects. In essence, SF 
within a cloud does not depend on how the molecular gas reaches the critical density, only that 
the density is reached before the SF clumps are disrupted by outside forces ( Elmegre"enl|201 1 ) . Our 
analysis of SFR within G38.9-0.4 appears consistent with this hypothesis. 
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6. Conclusions 

We have performed a multiwavelength analysis on the G38.9-0.4 SFR, which hes along a 
relatively unconfused line of sight at a distance of 2.7 kpc. The complex hosts a series of IR-bright 
bubbles, a web of IRDCs, and a large molecular cloud. 

1. The maximum gas mass surface density within the region is 0.2 g cm^-^. This suggests that 
current SF will is limited to low- and intermediate-mass stars, consistent with observations 
of the YSOs. 

2. The N 74 and N 75 bubbles are each powered by a single 09.5V star surrounded by a compact 
cluster. There is a similar cluster in MCM 16, but likely anchored by an early-B star. 

3. We detect a total of 162 YSOs with 54 stage I, 70 stage II, 6 stage III, and 32 ambiguous 
classifications. The high ratio of stage I to stage II YSOs (0.76) is consistent with the idea 
that SF in this region has only recently begun. 

4. The calculated YSO mass function slope a = —3.8 is significantly steeper than Salpeter but 
consistent with other YMF studies. 

5. We find only a weak correlation between individual YSO mass and the local gas mass surface 
density. 

6. There is, however, a strong correlation between areal YSO mass surface density and gas mass 
surface density with a tight power-law slope of 1.3. 

Classical signposts of triggering (pillars, YSOs coincident with bubble rims &: shells, YSO 
age gradients within bubbles) appear to be absent in this region. We conclude that, at present, 
traditional mechanisms (CC, RDI) proposed as triggers for SF have likely not occurred. This 
may be the case because either the single 09.5V stars provide insufficient mechanical/luminous 
energy to initiate triggering, or because the bubbles are too young for triggering to have begun, 
but we cannot discriminate between these possibilities at present. We used the bubbles to divide 
G38.9-0.4 into possible causal (triggered) and non-causal (non-triggered) regions to investigate the 
SF properties separately. We find no evidence to suggest that stellar feedback has enhanced or 
supressed the local SFE within the bubbles relative to other regions of the molecular cloud. This 
suggests that the SF rate and SFE are sensitive mainly to the gas density, regardless how the 
gas is collected. Detailed studies of the molecular gas around a larger sample of H II regions 
should be able to quantify the importance of feedback-triggered SF by determining the fraction of 
high-density, feedback-compressed gas surrounding H II regions relative to that already present in 
quiescent molecular clouds. Such studies would yield a quantitative measure of the fraction of all 
SF possibly influenced by stellar feedback. 



23 



MJA is supported by NASA GSRP fellowship NNX 10-AMlOH. HAK and CRK are supported 
by NASA through grant ADAP-NNX10AD55G. The authors would hke to thank Tom Robitaille 
for help with the YSO fitter and Matt Povich for his guidance in classifying YSOs. Thanks to 
Ed Churchwell and Barb Whitney for many helpful comments. The manuscript was significantly 
improved by the comments of an expert referee. 



REFERENCES 

Andre, P., Ward-Thompson, D., &: Barsony, M. 2000, Protostars and Planets IV, 59 

Allen, L. E., Calvet, N., D'Alessio, P., et al. 2004, ApJS, 154, 363 

Alexander, M. J., & Kobulnicky, H. A. 2012, ApJ, 755, L30 

Arvidsson, K., Kerton, C. R., Alexander, M. J., Kobulnicky, H. A., & Uzpen, B. 2010, AJ, 140, 
462 

Ballesteros-Paredes, J., Hartmann, L. W., Vazquez-Semadeni, E., Heitsch, F., & Zamora-Aviles, 
M. A. 2011, MNRAS, 411, 65 

Balser, D. S., Rood, R. T., Bania, T. M., & Anderson, L. D. 2011, ApJ, 738, 27 

Battersby, C, Bally, J., Ginsburg, A., et al. 2011, A&A, 535, A128 

Beaumont, C. N., & Williams, J. P. 2010, ApJ, 709, 791 

Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2003, PASP, 115, 953 

Bertoldi, F. 1989, ApJ, 346, 735 

Billot, N., Noriega-Crespo, A., Carey, S., et al. 2010, ApJ, 712, 797 

Bonatto, C, & Bica, E. 2009, MNRAS, 392, 483 

Cardelh, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 

Carey, S. J., Noriega-Crespo, A., Mizuno, D. R., et al. 2009, PASP, 121, 76 

Churchwell, E., Povich, M. S., Allen, D., et al. 2006, ApJ, 649, 759 

Clemens, D. P. 1985, ApJ, 295, 422 

Dale, J. E., Ercolano, B., & Bonnell, I. A. 2012, MNRAS, 427, 2852 

Deharveng, L., Zavagno, A., & Caplan, J. 2005, A&A, 433, 565 

Deharveng, L., Lefloch, B., Kurtz, S., et al. 2008, A&A, 482, 585 



-24- 

Deharveng, L., Zavagno, A., Anderson, L. D., et al. 2012, A&A, 546, A74 

Elmegreen, B. G., & Lada, C. J. 1977, ApJ, 214, 725 

Elmegreen, B. G. 2011, EAS Publications Series, 51, 45 

Fukuda, N., & Hanawa, T. 2000, ApJ, 533, 911 

Gritschneder, M., Burkert, A., Naab, T., & Walch, S. 2010, ApJ, 723, 971 

Groenewegen, M. A. T. 2006, A&A, 448, 181 

Gutermuth, R. A., Myers, P. C., Megeath, S. T., et al. 2008, ApJ, 674, 336 

Gutermuth, R. A., Megeath, S. T., Myers, P. C., et al. 2009, ApJS, 184, 18 

Gutermuth, R. A., Pipher, J. L., Megeath, S. T., et al. 2011, ApJ, 739, 84 

Herwig, F. 2005, ARA&A, 43, 435 

Indebetouw, R., Mathis, J. S., Babler, B. L., et al. 2005, ApJ, 619, 931 

Indebetouw, R., Robitaille, T. P., Whitney, B. A., et al. 2007, ApJ, 666, 321 

Jackson, J. M., et al. 2006, ApJS, 163, 145 

Kang, M., Bieging, J. H., Kulesa, C. A., & Lee, Y. 2009, ApJ, 701, 454 

Kendrew, S., Simpson, R., Bressert, E., et al. 2012, ApJ, 755, 71 

Kennicutt, R. C., Jr. 1998, ApJ, 498, 541 

Koenig, X. P., Allen, L. E., Gutermuth, R. A., et al. 2008, ApJ, 688, 1142 

Koenig, X. P., Leisawitz, D. T., Benford, D. J., et al. 2012, ApJ, 744, 130 

Krumholz, M. R., & McKee, C. F. 2008, Nature, 451, 1082 

Kurucz, R. L. 1993, Kurucz CD-ROM, Cambridge, MA: Smithsonian Astrophysical Observatory, 
— cl993, December 4, 1993 

Lawrence, A., Warren, S. J., Almaini, O., et al. 2007, MNRAS, 379, 1599 

Lefloch, B., & LazarefF, B. 1994, A&A, 289, 559 

Lucas, P. W., Hoare, M. G., Longmore, A., et al. 2008, MNRAS, 391, 136 

Maia, F. F. S., Corradi, W. J. B., & Santos, J. F. C, Jr. 2010, MNRAS, 407, 1875 

Marcolino, W. L. F., Bouret, J.-C, Martins, F., et al. 2009, A&A, 498, 837 



-25- 

Marigo, P., Girardi, L., Bressan, A., et al. 2008, A&A, 482, 883 

Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 

Markwardt, C. B. 2009, Astronomical Data Analysis Software and Systems XVIII, 411, 251 

Matsakis, D. N., Evans, N. J., II, Sato, T., & Zuckerman, B. 1976, AJ, 81, 172 

Mercer, E. P., Clemens, D. P., Meade, M. R., et al. 2005, ApJ, 635, 560 

Molinari, S., Swinyard, B., Bally, J., et al. 2010a, PASP, 122, 314 

Molinari, S., Swinyard, B., Bally, J., et al. 2010b, A&A, 518, LlOO 

Ossenkopf, V., k Henning, T. 1994, A&A, 291, 943 

Povich, M. S., Churchwell, E., Bieging, J. H., et al. 2009, ApJ, 696, 1278 

Povich, M. S., & Whitney, B. A. 2010, ApJ, 714, L285 

Povich, M. S., Smith, N., Majewski, S. R., et al. 2011, ApJS, 194, 14 

Rathborne, J. M., Jackson, J. M., & Simon, R. 2006, ApJ, 641, 389 

Rathborne, J. M., Jackson, J. M., Chambers, E. T., et al. 2010, ApJ, 715, 310 

Reach, W. T., Rho, J., Young, E., et al. 2004, ApJS, 154, 385 

Rho, J., Reach, W. T., Lefloch, B., & Fazio, G. G. 2006, ApJ, 643, 965 

Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 2006, ApJS, 167, 
256 

Robitaille, T. P., Whitney, B. A., Indebetouw, R., & Wood, K. 2007, ApJS, 169, 328 

Robitaille, T. P., Meade, M. R., Babler, B. L., et al. 2008, AJ, 136, 2413 

Roman-Duval, J., Jackson, J. M., Heyer, M., Rathborne, J., & Simon, R. 2010, ApJ, 723, 492 

Salpeter, E. E. 1955, ApJ, 121, 161 

Samal, M. R., Pandey, A. K., Ojha, D. K., et al. 2012, ApJ, 755, 20 

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 

Simpson, R. J., Povich, M. S., Kendrew, S., et al. 2012, MNRAS, 424, 2442 

Simon, R., Jackson, J. M., Clemens, D. P., Bania, T. M., & Heyer, M. H. 2001, ApJ, 551, 747 

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 



-26- 

Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 319, 730 

Stil, J. M., Taylor, A. R., Dickey, J. M., et al. 2006, AJ, 132, 1158 

Tafalla, M., Myers, P. C, Caselli, P., & Walmsley, C. M. 2004, A&A, 416, 191 

Testi, L., Palla, F., Prusti, T., Natta, A., k Maltagliati, S. 1997, A&A, 320, 159 

Testi, L., Palla, F., & Natta, A. 1999, A&A, 342, 515 

Thompson, M. A., Urquhart, J. S., Moore, T. J. T., & Morgan, L. K. 2012, MNRAS, 421, 408 

Urquhart, J. S., Busfield, A. L., Hoare, M. G., et al. 2008, A&A, 487, 253 

Watson, C, Povicli, M. S., Churchwell, E. B., et al. 2008, ApJ, 681, 1341 

Weaver, R., McCray, R., Castor, J., Shapiro, P., & Moore, R. 1977, ApJ, 218, 377 

Whitney, B. A., Wood,K., Bjorkman, J. E., & Cohen, M. 2003, ApJ, 598, 1079 

Whitney, B. A., Sewilo, M., Indebetouw, R., et al. 2008, AJ, 136, 18 



This preprint was prepared with the AAS lATgX macros v5.2. 



27 



Table 1. G38.9-0.4 Region Parameters 
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Fig. 1. — Three-color image of the G38.9-0.4 region. Blue is UKIDSS K, green is IRAC 8 /im, and 
red is MIPS 24 /im. The specific regions discussed in the paper are labeled. Green circles indicate 
the areas used for gas mass and radio continuum calculations, while magenta circles indicate the 
apertures used for cluster identification. The upper-most and lower-most white crosses mark the 



center location of the candidate cluster in (Mercer et al. 2005), while the inner crosses mark the 



location of the IMSFR IRAS source and UCEC. 
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Fig. 2.— Tb versus Vlsr plot of ^^CO for the main regions of G38. 9-0.4. Data for MCM 15 is 
shown as the black short-dashed hne, N 75 the orange sohd hne, N 74 the blue dot-dashed line, 
the IRDC the green dotted line, and MCM 16 the magenta long-dashed line. Associated molecular 
gas lies in the range of 35-45 km s^^, while peaks at higher velocities are foreground/background 
clouds with no apparent spatial correlation. The individual peaks show a slight shift towards higher 
velocities from North (MCM 15) to South (IRDC). The tangent velocity for this longitude is 79 
km s~^. 
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Fig. 3. — Channel plot showing ^^CO contours overlaid on an IRAC 8 /^ni image. The contours 
run linearly from 0.5 to 10.5 K km s~^ in steps of 1.25 K km s~^, and are averaged every 2 km s~^ 
over the range 34-46 km s~^. The velocity transition is seen across the different regions. The 
maximum value in each panel is (from left-to-right, top-to-bottom) 0.4 K km s~^, 5.6 K km s~^, 
12.0 K km s-\ 10.7 K km s-\ 8.0 K km s"\ and 1.8 K km s"^ 
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Fig. 4. — Temperature map of the SFR. Green and magenta circles show the location of important 
features as in Figure [T| Black contours indicate temperatures from 15 to 45 K, in steps of 3 K. 
The highest temperatures are associated with the bubble rims, while the coolest temperatures fall 
over IRDCs. The black region near / = 38?88 at the bottom of the image is due to lack of data 
necessary to complete the temperature fitting. 
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Fig. 5. — Three-color image with Herschel 500 /xm in blue, IRAC 8 /um in green, and MIPS 24 /im 
/im in red. White contours outline the VGPS 1.4 GHz continuum. The minimum level is the mean 
background temperature plus two times the standard deviation, which is 13.15 K + 2(0.43) K = 
14.01 K. Subsequent contours are 4, 6, 8, 10 and 12a above the mean background temperature. 
The contours show peaks coincident with each bubble, and some emission over MCM 15 with no 
clear peak. 
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Fig. 6. — Expansion radius versus age plot. Solid horizontal lines are the bubble radii of 0.65 pc 
for N 75 and 1.11 pc for N 74, while the curved lines (dotted, dashed, and dot-dashed) are the H II 
region expansion radii for densities of 10^, 10^, and 10^ cm~^. 
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Fig. 7. — Four-panel figure of YSO SEDs. The upper-left plot shows the SED for Y133 (stage I), 
the upper-right shows Y112 (stage II), the lower- left shows Y157 (stage III), and the lower-right 
shows an YSO (ambiguous). The solid black curve is the best-fit YSO model, the grey curves show 
the range of well-fit models, and the long-dashed curve is the SED for the underlying star. 



46 






CD 






12 

14 
3.0 

2.5 

2.0 

1.5 

1.0 

0.5 

0.0 
-0.5 




AA <^| 



-H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 h- 



+ + 






t^*.X^^ .^ - 





++ J 

+ 



H 1 1 1 1 1 1 h- 



+ 1 


All 


Dill 


O amb 



+ + 



+ 



JL A 



A^ct+5 + 

A ^ ^^ \ 6 ^ 



+^ 



DD 



e 



□ 



R ^ 



» 




^2.5 

2.0 



in 

.5 ^ 



1.0 ^ 

ID 

0.5 :! 



^0.0 
-0.5 



-0.5 0.0 0.5 1.0 1.5 

[5.8] - [8.0] 



2 4 6 

[8.0] - [24] 



Fig. 8. — Four-panel CCD and CMD plots using different color schemes. Stars lacking IR excess 
are shown as light gray dots. YSOs are labeled as plusses, triangles, squares, and circles for stage 
I, II, III, and ambiguous YSOs, respectively. Arrows indicate reddening vectors with Ay = 15 
magnitudes. 
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Fig. 9. — Three-color image with UKIDSS K, IRAC 8 //m, and MIPS 24 ^um in blue, green, and 
red. Over plotted are YSOs with stage I, II, III, and ambiguous labeled as red plusses, green 
triangles, blue squares, and cyan circles. White contours delineate the gas surface density levels of 
0.01(1.75^), equalling 0.01, 0.0175, 0.0306, 0.0536, 0.0938, and 0.1641 g cm'^. The YSOs primarily 
cluster around the bubbles and IRDCs. 
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Fig. 10. — Left: Three-color image of MCM 15. Blue, green, and red represent UKIDSS J, H, 
and K, respectively. The magenta circle indicates the region used for cluster identification. Red 
plusses, green triangles, and cyan circles represent stage I, II, and ambiguous YSOs. The contours 
outline MIPS 24 //m emission and are spaced as: d = 20(1.5*), or 20, 30, 45, 67.5, 101.3, 151.9, 
227.8, 341.7, 512.6, 768.9 MJy sr'^ Right: CMDs (top panel) and CCDs (bottom panel) for the 
MCM 15 target area (left), field area (center), and after field star substraction (right). The solid 
curve shows a ZAMS, and the dotted curve shows a 1 Myr PMS isochrones, both reddened by Ay 
= 12 and placed at 2.7 kpc. The asterisk and filled diamond indicate the location of a 09.5 and 
an AOV on the ZAMS. 
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Fig. 11. — Same as Figure 10, but for the N 75 bubble. Note the isochrone has been reddened by 
Av = 8. 
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Fig. 12. — TcQ — R relation for the candidate ionizing star in N 75. Filled black dots show the range 
Kurucz of stellar model that are well-fit to the star 2MASS 19034727+0509409. The gray hnes 
indicate the Teg — R for observed O stars, and the gray asterisks mark the location of 09.5 spectral 
types for the different luminosity classes. The data are consistent with 2MASS 19034727+0509409 
being a late-0 type star. 
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Fig. 13. — Same as Figure 10 but for the N 74 cluster. 
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Fig. 14. — Tgfj — R relation for the candidate ionizing star in N 74, with colors and symbols as in 



Figure [12] The stellar SED fits to the bright star, 2MASS 19035684+0506370, are consistent with 
a late-0 spectral type. 
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Fig. 15. — Same as Figure 10 but for the MCM 16 cluster. 
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Fig. 16. — Plot of the YSO mass function as a cumulative histogram. The number of YSOs, Ny^o, 
is shown as the solid black histogram, the best fit slope for YSOs > 2.5 Mq is given by the solid 
grey line, and a scaled Salpeter IMF is plotted as the dashed grey line. The completeness clearly 
begins to drop between 2.0 - 2.5 Mq and becomes quite significant below ~ 2 Mq. 
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Fig. 17. — Plot of YSO mass versus 500 iJ,ioa mass surface density. Stage I YSOs are plus symbols 
(upper left), stage II are triangles (upper right), ambiguous sources are circles (lower left), and 
stage III are squares (lower right). Large symbols mark the arithmetic average YSO mass and gas 
mass surface density for each stage. 



56 



o 

CL 



>- 
W 



^500 (g cm ) 
0.0002 0.002 0.02 

100.0 

0.0 



1.0 



0.2 



E 

0.002 



500 (g cm ) 



0.02 



0.2 



u 

CL 



>- 



)8:A 



10.0 



.0 



0.1 



o 



■'YSO 



A 




■^YSO 



0.01E 



1.38 ^Jy 



■'YSO 






a ° 


All O 


1 + 


II A 


amb O 



10 100 1000 

^500 (Mg Pc"^) 



10 100 1000 

^500 (Mg Pc"^) 



Fig. 18. — Plot of YSO surface density versus 500 ^m mass surface density. The upper left shows 
the surface density of all YSOs, the upper right shows only stage I, lower left shows stage II and 
lower right shows ambiguous objects. Solid lines indicate power law fits to the data given by the 
equation in the upper left-hand corner of each panel. 



